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L  INTRODUCTION 


Electrically  and/or  physically  small  antennas  are  often  required  by  missile  systems 
because  of  limited  space  requirements  or  reduction  in  radar  cross  section.  For  some 
systems,  even  at  UHF  and  VHF  frequencies,  the  antenna  can  be  isolated  as  the  single  most 
heavy  ana  bulky  component  (Reference  1).  Reduction  in  size  presents  problems  for  the 
system  designer  due  to  the  performance  penalties  in  bandwidth  and  efficiency  that  arise. 
These  performance  penalties  must  be  absorbed  into  the  overall  system  performance 
resulting  in  poor  reception  in  low  signal  regions.  Matching  techniques  can  be  used  to  boost 
antenna  performance  but  a  reduction  in  total  efficiency  is  still  incurred. 

Several  methods  to  counteract  the  above  problems  have  been  considered.  For  small 
antennas  with  radiation  resistances  lower  than  the  ohmic  resistances  of  their  radiating 
elements  and  which  include  feed  and  matching  networks,  higher  efficiencies  have  been 
measured  when  all  components  are  composed  of  high  temperature  superconducting  (HTS) 
materials  (References  2  through  5).  However,  this  replacement  of  normal  conductors  with 
HTS  components  can  only  partially  compensate  for  the  drop  in  efficiency  upon  size 
reduction.  Thus  it  is  advantageous  to  obtain  as  large  an  efficiency  as  possible  from  a  given 
rntenna  by  ether  means  first. 

One  such  method  has  been  investigated  for  some  years  by  D.  M.  Grimes  (References  6 
through  11).  This  method  consists  of  using  judicious  combinations  of  electric  and 
magnetic  multipole  sources  to  attempt  to  exceed  the  small  antenna  limitations  derived  by 
Chu  (Refer  :nce  12),  Harrington  (References  13  and  14),  and  Wheeler  (References  15  and 
16).  Since  most  of  Grimes'  analyses  on  this  subject  are  very  general,  in  this  report  we 
investigate  the  most  simple  case  of  the  above,  i.e..  an  electrically  small  colocated  array 
composed  of  an  electric  dipole  and  a  magnetic  loop  (modelled  as  an  orthogonal  equivalent 
magnetic  dipole).  We  have  found  that,  even  for  this  simple  case,  the  individual  electric  and 
magnetic  field  components  combine  in  such  a  way  that  the  interaction  term  involving  the 
currents  on  both  array  elements  disappears  from  the  complex  power.  This  allows  the 
reactive  part  of  the  complex  power  to  be  set  to  zero  when  appropriate  amplitude  and  phase 
ratios  between  the  dipole  current  and  the  loop  current  are  used.  Inputting  the  derived 
currents  into  the  numerical  code,  NEC3D,  we  find  that  the  radiated  power  is  significantly 
increased  over  the  corresponding  power  resulting  from  an  electric  dipole  or  a  magnetic  loop 
alone.  We  have  determined  that  die  loop  current  must  be  driven  in  phase  quadrature  to  the 
dipole  current  in  order  for  this  increased  radiated  power  effect  to  be  seen.  If  the  currents 
are  simply  in  phase,  the  array  radiates  exactly  as  if  each  element  stood  alone.  Similarly  the 
directivity  of  the  array  is  the  same  as  that  of  a  dipole  alone  or  a  loop  alone  if  the  respective 
currents  are  in  phase,  but  the  directivity  of  the  array  is  doubled  if  the  respective  currents  are 
in  phase  quadrature. 

In  Section  33,  we  determine  the  total  field  components  of  the  electric  dipole,  magnetic 
loop  array  in  simplified  form,  i.e.,  without  resorting  to  Hankel  and  Legendre  functions. 
In  Section  III  we  compute  the  radial  part  of  the  complex  Poynting  vector,  the  radial 
complex  power,  and  the  current  amplitude  ratio  which  gives  zero  reactive  power.  In 
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Section  IV  we  derive  the  radiation  intensity  and  directivity  of  the  array  for  two  cases:  (1) 
when  the  dipole  and  loop  currents  are  in  phase  and  (2)  when  they  are  in  phase  quadrature. 
The  NEC3D  numerical  results  are  presented  in  Section  V,  a  general  discussion  on  antenna 
element  coupling  is  presented  in  Section  VI,  and  conclusions  appear  in  Section  VII. 


II.  THEORY 

Consider  the  geometry  of  Figure  1.  We  assume  that  the  antenna  consists  of  two 
elements,  an  electric  dipole  oriented  along  the  z  axis  and  a  magnetic  loop  oriented  in  the 
same  plane,  the  yz  plane.  The  loop  and  dipole  are  colocated  and  nontouching.  The  dipole 

runs  from  -  U2  to  ill  along  z  with  a  total  length  i  and  a  current  Id.  The  loop  has  a  radius, 
"a,"  and  a  current  1/  We  assume  that  both  elements  are  electrically  very  small  and  thus 

considered  to  be  infinitesimal  (<X/50)  for  analysis  purposes. 

z 

A 


^  y 


FIGURE  1.  Antenna  Geometry. 
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Since  the  dipole  is  oriented  along  the  £  axis,  its  electromagnetic  field  components  in 
spherical  coordinates  (for  a  constant  current)  are  (Reference  1 7): 


^  4  7tr  V  lAr/ 


£(i)  =  5M£2L£i1+-LV^ 

'  2wr  v  ikrj 


(la) 

(lb) 

(lc) 


_  irjklji  sin  $ 


4  nr 


i+J — 


1 


ikr  (kr) 


7 


(Id) 


where  k 


2n  _  [a 

=^,n  v? 


and  the  superscript  1  indicates  that  these  are  dipole  fields.  Now 


we  must  put  the  loop  into  the  same  coordinate  frame.  Standard  analyses  of  loop  antennas 
place  them  in  the  xy  plane  with  £  normal  to  the  loop.  Tn  our  case  the  loop  is  located  in  the 
yz  plane  with  £  as  its  normal.  We  invoke  the  usual  equivalency  between  an  infinitesimal 
loop  of  radius  a,  current  If,  and  an  infinitesimal  magnetic  dipole  of  length  current  Im 


which  (in  our  case)  runs  along  the  £  axis.  Thus  (Reference  17) 


Vm  =  i(na2)o)ki/f  =  i(mi2)ki]lf  . 
An  equivalent  magnetic  dipole  along  the  £  axis  has 


(2) 


A  =0 


(3) 


F(x'y’z)  =  J^! 


r  .  r  _r\  C 


-IkX 
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- df 
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where  A  and  F  are  the  usual  magnetic  and  electric  vector  potentials,  respectively.  Using 

Im-  Im^  where  Im  is  assumed  to  be  constant  and  assuming  that  R  «*  r,  the  electric  vector 
potential  becomes 


~p  -  % 

4  fir 

Using  the  transformation  between  Cartesian  and  spherical  coordinates, 


(4) 


(Fr) 

rsin$cos^  sin  8  sin  <f>  cos  8  s 

(px) 

Fe 

= 

cos  6  cos  <p  cos  8  sin  ip  -sin  8 

0 

(5) 

K  -sin  0  cos  <p  0  , 

,0; 

we  obtain 


_  *  /* 

Fr  =  -—^  sin  6  cos  <f>  e~ikr  (6a) 

4  7TT 

Fq  -  cos  8  cos  <p  e~ikr  (6b) 

A  nr 

F  ~e/ym  sin0e-.*r  (6c) 

Afir 

The  electric  field  components  of  the  loop  are  given  by 


or 


(7) 


(8a) 
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42)  =  -^a^-sin  J 1  + — \~ikr 

6  Aicr  Y{  ikr  J 


(8b) 


£(2)  _  IMmim cos  Q  cos  lj>l  1  -f  —If  >kr 
*  4  nr  \  ikr) 


(8c) 


where  the  superscript  2  indicates  that  these  are  loop  fields. 
Using 


H  -  ——VxE  , 
icon 


(9) 


die  magnetic  field  components  are 


h(2)  ^  sjn  0  cos  </  1  +  •— Vikr  . 

r,2zr *  ‘  V  /Ary' 


Hfial 


/*>\  —  »it/  / 

//<‘)=.^fffi£filcos0cos0 

e  r/47rr 


1  1 


1  +  — - 


(Ar  (kry 


-ikr 


(10b) 


isin0 

*  r,4;m  W 


1  + 


i  >  u 

ikr  {krf[ 


(10c) 


Thus  the  total  field  components  of  the  dipole  and  the  loop  should  be  simply  a  superposition 
of  the  field  components  from  each  element  alone.  Using 


E  =  E{l)+E(2)  ;  H  =  Hw+H{2)  , 


(11) 


the  total  field  components  are 
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£r  =  -5i 4cos  e{)  +  — 

r  2k?  I  ikr) 


(12a) 


Eq  -  lJ]Mdis  w  Q 
0  4  nr 


1  +  - ^TV1'*r  +  i^o-sin  *(  1  +  — V*r 
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Equations  12  are  a  suf-erposibon  of  a  TM  (to  ?)  mode  from  the  dipole  and  a  TE  (to  i?)  mode 
from  the  loop.  The  field  components  from  the  dipole  have  no  <}>  dependence,  while  those 
from  the  loop  have  both  0  and  <J>  dependence  in  genera)  (although  and  are  not 

functions  of  0)  Of  the  total  field  components  given  by  Equations  12.  none  are  zero,  but 
only  two,  Ee  and  H^,  contain  contributions  from  both  the  dipole  and  the  loop. 


III.  THE  RADIAL  PART  OF  THE  COMPLEX  POINTING  VECTOR 

AND  COMPLEX  POWER 


The  usual  complex  Poynnng  vector  is  defined  to  be  (References  17  and  18) 


(13) 
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or  in  component  form 

at, * 

SCu=\[E0Hr-ErHl)  , 
SCt=±[ErH6-EeHr)  . 


(14a) 


(14b) 


(14c) 


For  now,  we  will  consider  only  the  radial  pan  of  the  complex  Poynting  vector  given  by 
Equation  14a.  Using  Equations  12  and  after  some  algebra,  Equation  14a  becomes 


Sc  =  - 
c'  2 


* 

1  "  T 

(4^r)2 

(kry 

4lVi^)-!singsir,0 

(4  »r 

+  sintf  sin  » 

(4*r)2 


,  1  1 

1 - >r  + 


(*rf  (U?"J 

1  +  A 

(*r)2 


2  |2 

+  iM4_[sin2  0  +  cos2  e  cos2  0] 

1  +  — U- 

r?(4/tr)2  1  J 

{krf  J 

(15) 


To  obtain  complex  power  travelling  in  the  radial  direction,  we  integrate  the  radial  part  of  the 
complex  Poynting  vector  over  a  sphere  of  radius  r.  So 
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r  =  §  Sc  d~*s  (16) 

In  k 

=  j  J  [V  +  Sc>  +  V]  [r2  sin  ed6d(pr] 

0  0 

or 

In  n 

P  =  r‘  J  J  SCr  sin  0  d$d$ 

0  0 

with  Scr  given  by  Equation  15.  Considering  the  four  terms  in  Equation  15  and  since  we 
only  have  to  integrate  over  the  angular  dependence,  the  first  term  in  Equation  15  leads  to  an 
integral  of  the  fonn 

2/r  n 

/,  =  j  J  sin3  d  dOd<p  ,  (19) 

0  0 

the  second  and  third  terms  in  Equation  15  lead  to 
2  n  x 

h  =  j  j  s'n2  #  sin  0  dfkUp  ,  (20) 

0  0 

and  the  fourth  term  in  Equation  15  leads  to 
2»r  n 

1$  =  J  J  (sin2  0  +  cos2  0cos?  0)si»  9d6d<p  .  (21) 

o  o 

Performing  the  integrations  in  Equations  19  through  21,  Ij  =  y,  I2  =  0, 13  =  Thus  the 
complex  radial  power  becomes 


(17) 


(18) 
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P  = 


Vk2\ldt\2 

\  1  , 

1  +  -U-1 

(4  n)2 

(*r)5_ 

r?(4/i)2 

L  (Ar)3 J 

(22) 


Because  the  second  integration,  given  by  Equation  20  vanishes,  the  interaction  terms 
between  the  dipole  and  the  loop  disappear  from  the  radial  complex  power.  These  cross 
terms  are  still  present  in  the  complex  Poynting  vector. 

Splitting  Equation  22  into  real  and  imaginary'  parts, 


Re(P)  = 


(23a) 


and 


Im  (/>)  = 


( 


1“  \ 


1 

(*r7 


(23b) 


If  we  interpret  Equation  23a  as  the  radial  radiated  power  and  Equation  23b  as  the  radial 
reactive  or  stored  power,  then 


rrad  = 


r{n 


3A 


m^m 


\2\ 


(24a) 


and 


(24b) 


where  Prad  is  the  time-averaged  radiated  power,  cu  is  the  frequency,  and  Wm  and  We  are  the 
time-averaged  magnetic  and  electric  radial  energies,  respectively. 

Due  to  the  disappearance  of  the  interaction  terms  between  the  loop  and  the  dipole  in 
the  complex  radial  power,  we  can  set  the  reactive  or  stored  power  to  zero  by  insisting  that 
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/  / 
X 


T 


(25) 


Using  Equation  2,  the  proper  current  amplitude  ratio  for  zero  reactive  power  is 


However,  Equation  26  does  not  tel!  us  about  the  relative  phases  of  the  dipole  and  loop 
currents.  We  will  discuss  this  issue  in  Section  IV. 


IV.  RADIATION  INTENSITY,  DIRECTIVITY,  AND 
RELATIVE  CURRENT  PHASES 


From  the  total  field  components  in  Equations  12,  in  the  far-field  these  components 
become 


Er~Hr*>  0  , 


(27a) 


-Her 


Eq  ~  — — - ( riljt  sin  6  +  iAIdf  sin  <p)  , 

4  tit 


(27b) 


9  4/rr  v 


(27c) 


with 


He  =  -r)E<)  ;  H<p  =  T)Eq 


(28) 


and  A  =  7ta2,  the  loop  area. 

The  radiation  intensities,  Uq  and  U$,  are 
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Ue  = 


rl\Ee 

2  r. 


u*~ 


lT\ 


(29) 


and  wc  will  consider  two  separate  cases  based  on  the  phase  relationship  between  the  loop 
and  dipole  currents. 

Case  1.  Let  I<i  be  pure  real  and  be  pure  real.  Then  since  they  are  in  phase 


U9  =  * 
0  8 


sin2  #  + 


(2; 


i  uJ 


\Uj 


sin2  <p 


(30a) 


and 


Ua  ~ 


W\T  I 

Jd) 


9  9 

cos  6  cos  <P 


(30b) 


and  the  total  intensity  is 


d  cos2 


(31) 


Since  the  radiated  power  is  given  by  Equation  24a,  wc  can  obtain  both  the  directive  gain 
from 


D  -1^ 

*  prad 


and  the  directivity  from 


D  - 

°  Prod 


(32) 


(33) 
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The  directive  gam  is  plotted  for  4>  =  ^ J  *n  Figure  2  with  ^ - J  as  a  parameter  for  I;  and  I<j  in 

phase.  The  maximum  directivity  of  this  antenna  (when  the  currents  are  in  phase)  is  3/2 
(1.76  dB),  exactly  the  same  value  as  the  directivity  of  the  small  dipole  alone,  or  the  small 
loop  alone  (Reference  17). 


IjL  and  I,j  are  in  phase 


FIGURE  2.  Directive  Gain  When  Array  Element  Currents  are  in  Phase. 

Case  2.  Let  1^  be  real  and  L  =  -il/  where  1/  is  real.  The  two  currents  are  now  90  degrees 
out  of  phase,  or  in  quadrature.  Now 


-ikr 


Eq  ~  - (Ijf  sin  6  +  Akl(  sin  0) 


(34a) 


and 


„  ik2A-nl'fe~ikT 
- - 


Anr 


cos  0  COS  <P 


(34b) 


Thus  the  radiation  intensities  for  this  case  are 
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H 

\!dj 


sin  6  sin  0 


.  4  a 
+Ait  -- 

\  A 


f©  f 


2  2 

sin  0+cos 


d  cos" 


0) 


(35a) 


and 


(35b) 


Immediately  we  see  that  Equation  35a  contains  extra  terms  when  compared  with  Equation 
30a.  U(j>  as  given  by  Equations  30b  and  35b  is  the  same  for  both  cases.  When  <j>  =  0, 

Equations  30a  and  35a  are  identical  for  the  same  | j  When  $  =  Equation  35a  is  larger 

tlmn  Equation  30a  (see  Figure  3).  Where  the  directive  gain  is  computed  using  Equation  35, 
the  amplitude  as  well  as  the  shape  of  the  curves  are  quite  different  from  the  in  phase  case. 
_  !/ 

The  maximum  directivity  is  3(4.77  dB)  in  this  case  when  l-H  is  given  by  Equation  26. 

Thus  the  antenna  directivity  for  currents  in  quadrature  is  twice  as  much  as  the  directivity 
resulting  when  the  currents  are  in  phase. 
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Ij  and  1^  are  in  quadrature 


V,  MUTUAL  AND  DRIVING  POINT  IMPEDANCES  USING  NEC3D 


This  section  describes  the  approach  used  in  modelling  the  colocated  magnetic  loop 
and  electric  dipole  array  shown  in  Figure  1  (also  sec  Reference  19).  The  modelling 
approach  can  be  summarized  by  the  following  steps: 


( 1 )  Find  the  array  Y-matrix  using  Method  of  Moments  (MOM) 

(2)  Find  the  array  Z-matrk  (by  inversion  of  Y-matrix) 

(3)  Determine  desired  current  excitations  for  both  loop  and  dipole 

(4)  Determine  active  port  impedances  for  both  loop  and  dipole  from: 


za, 


=  h\  f ■  +  zi2-f+-+zu+-+zin  -f- 


h 


h 


a 


(36) 


where 

Zai  is  the  active  impedance  at  the  ith  port 
Zij  is  the  ijth  element  of  the  Z-matrix 
I,  is  the  current  at  the  ith  port 

The  remainder  of  this  section  will  discuss  each  of  the  above  steps  in  more  detail  and  present 
results  for  the  colocated  array. 
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FINDING  THE  Y-MATRIX  USING  MOM 

In  the  MOM  approach  the  antenna  array  geometry  is  first  defined  by  wire  segments  as 
shown  in  Figure  4.  The  electrically  small  loop  to  be  analyzed  has  been  assumed  to  be 

square  with  side  length  of  0.025  X.  Each  side  of  the  loop  is  subdivided  into  five  segments. 
The  loop  is  centered  about  the  origin  and  is  located  in  the  yz  plane.  The  dipole  is  also 

centered  about  the  origin  and  is  located  along  the  z-axis  with  total  length  of  0.020  X.  The 
dipole  is  also  divided  into  five  segments.  For  this  example,  the  radius  of  the  wire  was 

chosen  as  0.001  X.  The  total  length  of  the  dipole  was  adjusted  to  be  80%  of  the  loop  side 
length  to  avoid  close  proximity  of  the  ends  of  the  dipole  with  segments  in  the  loop.  The 
effects  of  varying  the  relative  dipole  length  have  not  been  investigated  but  changes  are  not 
expected  to  change  the  basic  operating  principles  of  the  colocated  array. 


z 


FIGURE  4.  Colocated  Antenna  Geometry  for  Method  of  Moments  Analysis. 


In  MOM  the  array  is  excited  by  inserting  delta  gap  voltage  generators  of  specified 
magnitude  and  phase  at  the  center  of  a  wire  segment.  The  cunrents  arc  then  determined  in 
each  of  the  segments  according  to  the  MOM  procedure.  As  wili  be  shown  below  the  MOM 
approach  lends  itself  to  direct  computation  of  the  two  port  Y  matrix  of  the  array.  The  Y 
matrix  of  a  two  port  is  described  in  terms  of  the  terminal  voltages  and  currents  by: 


15 


NAWCWPN'S  TP  8212 


11  =  Yn  V1  +  Y12V2 

12  =  Y2t  Vi  +  V22  V2  .  (37) 


Notice  that  if  Vj  is  set  to  1  volt,  and  V2  is  set  to  zero:  Y 1 1  is  identically  equal  to  “he  current 
li  and  Y21  is  equal  to  I2.  Similarly,  if  Vj  is  set  to  zero  and  V2  is  set  to  1  volt,  then  Y  12  is 
equal  to  the  current !)  and  Y22  is  equal  to  I2.  In  this  fashion  the  Y  matrix  can  be  directiv 
computed  from  the  MOM  detived  currents  since  the  delta  gap  voltage  can  be  set  to  1  volt 
and  the  voltage  set  to  zero  by  replacing  the  delta  gap  generator  by  a  wire  segment  (i.e. 
removing  the  delta  gap  generator).  The  Y  matrix  computed  for  the  example  in  Figure  4 
(with  lossless  conductors;  at  a  frequcnc y  of  500  MHz  is: 


f  1.1293  - 10-9 

1^7.7687  •  10~9 


3.41 1 1  •  KF4/ 
5.5983  10  _5i 


7.4836  - 10-9 
1.4187- HT6 


■  5.8541 -HTV 
0.01039/  J 


(38) 


The  above  values  were  computed  using  Numerical  Electromagnetics  (..'ode  (NEC), 
version  NEC3D  run  on  a  Macintosh  Ilsi.  NEC  was  developed  at  the  !_awrence  Livermore 
Laboratory  under  the  sponsorship  cf  tne  Naval  Ocean  Systems  Center  and  the  Ah  Force 
Weapons  Laboratory- 


ARRAY  Z  MATRIX 

The  Y  matrix  can  be  inverted  to  give  the  Z  matrix 

f  0.00934- 2.9289 -103/  9.13542  •  10 '5 +  16.50405/^1  ...... 

1,-8.5203.?  •  1(T5  4- 15.78289/  0.01314  +  96.16672/  J 

Zn  corresponds  to  the  self  impedance  of  the  dipole  and  Z22  is  the  self  impedance  of 
the  loop.  Since  in  this  example  conductors  are  assumed  lossless,  Rj  1  and  R22  represent 
the  radiation  resistances  of  the  dipole  and  loop,  respectively,  and  since  both  loop  and  dipole 
are  electrically  small,  the  radiation  resistance  of  each  isolated  antenna  is  correspondingly 
very  small.  One  of  the  prime  motivations  for  this  hybrid  mode  array  is  to  significantly 
increase  the  radiation  resistance  of  the  electrically  small  antennas  without  increasing  their 
physical  size. 


COMPUTING  CURRENT  EXCITATIONS  OF  HYBRID  MODE  ANTENNA 

Using  Equation  26  from  Section  III  that  specifies  the  proper  current  ratio  for  zero 
reactive  power  and  rewriting  the  denominator  in  terms  of  the  area  of  the  loop.  A,  gives 
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Id 


(40) 


For  a  square  loop  with  side  length  S  the  proper  ratio  is 


h 

Id 


(41) 


for  (/X  =  0.020,  S/1.  ■-  0.025,  and  I^/Id  =  i  5.093.  Thus,  the  current  in  the  loop  is 
approximately  five  times  that  in  the  dipole,  and  the  currents  are  in  quadrature. 


ACTIVE  IMPEDANCES  OF  HYBRID  MODE  ANTENNA 

The  active  impedances  are  computed  from  Equation  36  and  listed  in  Table  1  for  the 
relati  ve  current  excitations  shown. 


TABLE  1. 


dipole 

I|9?D 

^diDOle 

ZaiQOD 

1  Z0° 

5.093 

0.0098  -i2844.S 

0.01313 +  i99.266 

1  Z  90° 

5.093 

84.064  -i2928.9 

-3.086  +  i96.17 

1  Z  180° 

5.093 

0.0089  -  i301 3.0 

0.01316  +  193.068 

1  Z270p 

5.093 

-84.045  -  i2928.9 

3.112  +  i96. 167 

It  is  important  to  note  that  when  the  currents  are  in  quadrature  the  radiation  resistance 
of  both  the  loop  and  dipole  is  sigmficar  Jy  increased.  If  the  currents  are  not  in  quadrature 
the  radiation  resistance  is  approximately  that  of  the  isolated  element.  Furthermore,  with  the 
currents  in  quadrature  the  radiation  resistance  of  one  element  is  positive  and  the  radiation 
resistance  of  the  other  element  is  negative.  The  physical  interpretation  of  a  negative 
radiation  resistance  is  that  rather  than  radiating  power  (as  would  be  the  case  for  a  positive 
radiation  resistance),  an  antenna  with  a  negative  radiation  resistance  is  receiving  power  and 
delivering  it  to  the  load  connected  to  the  port.  Tne  significance  of  colocation  and 
quadrature  currents  in  the  large  increase  in  radiation  resistance  can  be  seen  in  the  expression 
for  the  active  impedances  of  a  two  port  colocated  antenna  array. 
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z&\  =zn  +  zi2 


ll 

h 


Zdi  —  Zj\ 


h-  +  Z22 


(42) 


From  Equation  39  the  mutual  coupling  (Z12  and  Z21)  for  colocated  antennas  is 
observed  to  have  a  large  reactive  component.  As  the  separation  between  antennas  is 
increased  the  magnitude  of  the  mutual  coupling  also  decreases.  The  large  reactive 
component  of  the  mutual  coupling  combined  with  the  currents  in  quadrature  results  in  a 
significant  real  resistive  component  in  the  active  impedance  compared  to  the  small  radiation 
resistance  due  to  the  self  impedance  of  electrically  small  antennas. 


SENSITIVITY  TO  MAGNITUDE  AND  PHASE  PERTURBATIONS 

Table  2  lists  computed  active  impedances  for  current  excitations  that  deviate  from  the 
specified  current  ratio  of  Equation  41. 


TABLE  2. 


mmmmsmm 

_ ^diPQle _ 

£aiooo 

!  P dipole  .  ^Iood 

1  Z  90° ,  6.000 

mi  “kn 

1  Z  90°  ,  5.093 

84.064  - 12928.9 

-3.086 +  i96. 17 

42.03  ,  -40.02 

1  Z  90°  ,  5.000 

82.531  -  i2928.9 

-3.144  -t- 196. 17 

41.27  ,  -39.29 

1  Z  90°  ,  4.000 

66.027  -  i2928.9 

-3.933  +  i96.17 

33.01  ,  -31.46 

;  Z  85° ,  5.093 

83.746  -i292 1.6 

-3.074  +  i96.44 

41.87  .  -39.87 

1  Z  45°  ,  5.093 

59.446  -  i2869.5 

-2.178  +  i98.36 

29.72  ,  -28.25 

1  Z  135°,  5.093 

59.446  -  i2988.4 

-2.178  +  i93.98 

29.72  ,  -28.24 

Thus,  it  appears  that  the  colocated  antenna  is  relatively  insensitive  to  small  perturbations  In 
the  magnitude  and  phase  of  the  excitation  currents.  The  powers  listed  for  the  dipole  and 
loop  are  the  radiated  power  from  the  dipole  (positive)  and  the  power  received  by  the  loop 
(negative),  respectively.  The  total  radiated  power  of  the  array  is  die  algebraic  sum.  As  die 
magnitude  of  the  current  in  the  loop  is  decreased  the  total  radiated  power  also  decreases. 
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VI.  GENERAL  DISCUSSION  ON  ANTENNA  ELEMENT  COUPLING 


In  antenna  arrays  the  relative  input  currents  to  the  individual  antenna  elements  are 
specified  as  to  their  relative  magnitudes  and  phases  in  order  to  obtain  the  desired  gain, 
pattern,  radiation  direction,  etc. 

With  close  packed  arrays,  (and  a  colocated  loop  and  dipole  is  certainly  an  example) 
where  the  coupling  between  the  elements  cannot  be  ignored,  a  serious  efficiency  problem 
can  arise. 

Consider  our  two-element  array  as  a  two  port  as  in  Figure  5a,  with  impedance  matrix 
(Z)  and  scattering  matrix  (S).  It  follows  that  the  active  input  impedances  looking  into  ports 
(or  antenna  elements)  one  and  two  are: 


Zi  -Zu  +  Z12  — 


(43a) 


Z2  -  Z22  4  Zl2 


Ll 


(43b) 


Since  the  Z,j  are  tixed  by  the  geometry  of  the  situation  and  I2  and  lj  arc  fixed  by  the 
antenna  performance  requirements,  Z\  and  Z2  are  fixed  impedances  with  specific  numerical 
values. 


The  active  impedances,  Zj  and  Z2,  represent  the  impedances  that  must  be  matched 
into  in  order  to  increase  the  system  efficiency  and  radiate  as  much  of  the  available  power  as 
possible.  The  matching  problem  reduces  to  Figure  5b  with  some  provisos  that  we  will  note 
as  we  proceed. 


t>2 


a2 


FIGURE  5(a).  Two-Element  Array  as  a  Two  Port. 


aA1 

bAl 


1  <SA)  2 


®1 


b1  •*--  u 


r. 


(SB)  1 


aBl 


bBl 


FIGURE  5(b).  Reduction  of  Two  Port  Matching 
Problem. 
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The  active  reflection  coefficients,  I'i  and  T2,  are  given  by 


r  .-hzlSL 

'  Zi+Zo 


i  - 1,2 


(44) 


while  (Sa)  and  (Sb)  represent  matching  networks  with  the  goal  of  reducing  the  reflected 
waves  bAi  and  bei  to  zero  while  simultaneously  not  introducing  loss  in  the  matching 
networks  themselves. 

For  a  lossless  matching  network,  the  condition  for  zero  reflection  is  given  by 


•Sa22  =  1'i  ;  5fl22=r2  *  <45) 

The  other  scattering  coefficients  of  the  matching  networks  can  then  be  found  by 
imposing  the  usual  unitary  condition 


(S)(sr)*  =((/)  (46) 

on  the  matching  networks. 

It  should  be  noted  that  merely  attaching  the  matching  networks  thus  designed  to  the 
antenna  array  input  ports  will  not  ensure  proper  performance  of  the  antenna  array.  The 
drive  levels  aAi  and  aei  into  the  matching  networks  must  be  of  the  proper  relative 
magnitude  and  phase  to  provide  the  correct  relative  currents  Ij  and  I2.  The  relations 
between  asi  and  the  SAij  and  current  Ij  are  readily  found  as  are  the  relations  between  am 
and  the  Sgij  and  current  I?. 

The  drive  levels  aAi  and  aBl  would  generally  be  supplied  from  an  uneven  power 
divider  designed  to  give  the  proper  complex  ratio  aAi /aBl  with  the  system  looking  like 
Figure  6. 
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FIGURE  6.  Two-Element  Array  With  Matching  Networks 


As  long  as  Equation  45  can  be  satisfied,  there  is  no  reason,  in  principal,  why  lossless 
matching  networks  and  power  dividers  cannot  be  designed  and  buiit  to  a  good  degree  of 
approximation.  If  the  coupling  between  antenna  elements  is  sufficiently  strong,  and  the 
currents,  Ij  and  I2,  have  sufficiently  different  magnitudes  and  phases,  the  real  pan  of  Z\  or 
Z2  (but  no;  both  due  to  the  conservation  of  energy)  can  be  negative.  When  that  is  the  case, 
the  magnitude  of  the  active  reflection  coefficient  is  greater  than  unity  and  cannot  be  matched 
using  a  lossless  network. 


For  the  electric  dipole  and  the  loop  with  Ij  and  lj  set  to  "zero  out”  the  imaginary  pan 
of  the  complex  Poynting  vector,  this  is  precisely  the  case  as  shown  in  Table  1.  The 
mechanism  by  which  this  occurs  is  readily  seen  by  rewriting  Equations  43a  and  43b: 


Zi  =  *j+  j  + 


~2-(/?i2cos  A<p-Xi2  sin  A$) 


+  1 


41 


+ 


(/?12  sin  A<p  +  X12  cos  A <j>) 


(47a) 


Z2  ~  ^2  ^2  =  ^ 22 


(/?! 2  cos  A<p  +  X12  sin  A^>) 


+  t 


^22  +  ^12  cos  A0-/?i2sin  A<p) 

h 


(47b) 
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where  A<j>  is  the  phase  shift  between  Ij  and  I2.  If  Il2l»llil.  for  example,  and  R12  is  small, 

and  A<J>  =  jt/2,  it  is  easily  seen  that  Ri  may  very  well  be  negative,  while  R2  is  positive. 
(R  j  1  and  R22  in  the  impedance  matrix  must  be  positive  without  the  presence  of  active 
devices.) 


If  a  negative  Rj  is  the  case  then  the  matching  problem  leads  to  Figure  7,  where  Fa  is 


chosen  to  give  the  proper  current  at  port  1  of  our  antenna  array. 

al 

1  (S)  2 

—  h2  U 

rA 

■  _ 

K 

M 

*2  -  - 

(%)  1 


*8 1 


FIGURE  7.  The  Matching  Problem  When  Rj  is  Negative. 


Figure  8  is  equivalent  to  Figure  7  and  maintaining  the  proper  value  of  Ij  must  be  done 
by  setting  Ta  equal  *0  the  reciprocal  of  fi. 


* —  "2 

— ►  h 

2  (SB)  1 

■* —  aBt 

r 

*2 

<2 

FIGURE  8.  Equivalent  System  to  Figure  7. 


=  i-6.)=|L(i-r,)=|L 

4)  4)  A) 


1- 


(48) 


Solving  (48)  for  I'a  gives 


(49) 
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(50a) 

where 

Z,  =-)/?i  |  +  iX,  .  (50b) 


Because  IRjl  is  a  positive  ohmic  resistance  and  the  current  Ii  flows  through  it,  there  is 
ohmic  loss  which  means  that  although  the  antenna  performance  in  terms  of  pattern  and 
radiation  direction  is  preserved,  the  efficiency  is  reduced.  In  the  case  of  electrically  small 
superconducting  antennas,  this  reduction  in  efficiency  can  be  drastic. 

However,  all  is  not  lost.  The  load  Za  can  be  replaced  with  a  lossless  two  port  as  in 
Figure  9.  The  lossless  matching  network  (Sg)  is  designed  for  matching  so  that 


SB22  =  r2 


FIGURE  9.  Replacing  the  Load  Za  With  a  Lossless  Two  Port. 


Per  Equation  45,  while  the  lossless  two  port  (Sa)  is  designed  to  produce  the  current  lj  with 


SA22  ~  TT  ’ 

1  1 


(51) 


the  output  bAi  is  fed,  with  the  proper  phase  shift,  to  the  input  of  a  power  combiner  to 
produce  the  desired  drive  level  agi- 

It  is  not  our  purpose  to  go  into  the  combiner  details  here  but  one  way  to  achieve  this 
end  is  to  use  a  directional  coupler  with  one  port  terminated  in  Zo,  with  the  coupling 
coefficients  set  to  give  zero  output  at  the  terminated  port  and  the  proper  agi  of  the  input  of 
the  matching  network. 

This  may  seem  to  be  a  somewhat  complex  process  but  if  it  is  compared  with  Figure 
6,  one  power  divider  and  two  two  port  matching  circuits  versus  one  power  combiner,  one 
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two  port  matching  circuit,  and  one  two  port  to  meet  the  condition  of  Equation  51,  the  two 
cases  are  seen  to  be  topologically  equivalent. 


VII.  CONCLUSIONS 


We  have  presented  a  detailed  electromagnetic  analysis  of  an  electrically  small 
colocated  electric  dipole  and  magnetic  loop  antenna  array.  This  antenna  is  the  simplest 
examole  of  the  Grimes  multipole  class  of  antenna  arrays.  We  have  determined  that  since 
the  interaction  term  between  the  two  elements  disappears  from  the  radial  complex  power, 
we  were  able  to  set  the  radial  reactance  to  zero  by  choosing  appropriate  current  magnitudes 
and  phases  on  the  array  elements.  By  driving  the  two  elements  in  quadrature,  we  saw  a 
much  increased  radiation  intensity  and  directivity  as  well  as  increased  radiated  power. 
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Appendix 


RELATIONSHIP  BETWEEN  THE  INSTANTANEOUS  POYNTING 
VECTOR  AND  THE  COMPLEX  POYNTING  VECTOR 


It  is  important  to  recall  that  the  complex  Poynting  vector  can  be  derived  directly  from 
Maxwell's  equations  (Reference  20).  Using  standard  phasor  notation  and  e‘wt  time 
dependence 


Vx  E  =  -i(Ofi  H 

(A-U) 

Vx  H  =  J  +  icoe  E 

(A- lb) 

with  the  complex  conjugate  equation  for  Equation  A- lb  as 


VxH*  =  J*-ico£F* 


(A-2) 


Dotting  P?*  into  Equation  A- la  and  ^  into  Equation  A-2,  we  obtain 


H*  VxE  =  -imH*  H 


(A-3) 


and 


E  VxH*  =  E  J*-io>£ E  E* 


(A  3b) 


Subtracting  Equation  A- 3b  from  Equation  A-3a,  we  have 
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->  A 


V^ExH*\  =  -ico\ 


-E  J*  . 


(A-4) 


Allowing  primed  variables  to  denote  real  parts  of  quantities  and  double  primes  to  denote 
imaginary  parts,  we  can  rewrite 


HH*-H'2  +H"2 


and 


(A-5a) 


E-  E*  =  E'2  +  E"2 


(A-5b) 


which  have  no  imaginary  parts. 

However,  E?  x  1?*  does  have  an  imaginary  part  Thus 


ExH*  =  \  E'x //'+  E" x H"  I  +  if  E'x H"~  E"xH ' ' 


Y 


(A-6) 


while 


£  ./  =E*  cr£’ 


(A-7) 


and  as  long  as  the  conductivity,  o,  is  real,  t  •  J  *  is  a  purely  real  quantity.  So  returning 
to  Equation  (A-4), 


/— >  — *  \ 

f  — >  - >  - >  - *  'N 

r— ►  — *  — *  —*  \ 

ExH*  -  V- 

E'xH’+  E"xH"Ui 

E'xH"-  E"  xW 

V  J 

A  y 

\  J. 

=  V 


f  —*  — i ►  )  f— *  ^ 

Re^Erf/  *J+iTm|^£r//  * 


(A-8) 
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Using  Equation  A-8  with  Equation  A-4, 


5V  Rl£^Tl[£'y*J 


(A-9a) 


and 


(  — »  — »  > 

(  — »  — >  -»  -t  \ 

H  H*  EE* 

i  ExH  * 

=  -0) 

fJ. — - - 

V  J 

2  2 

(A-9b) 


Integrating  Equations  A-9  over  a  sphere  of  volume  V  v/ith  boundary  S,  we  see  that 


EJ*dV 


(A- 10a) 


and 


V-Imf  =  -(yJ 

v'  '  '  v 


\ 


H  H*  E  E * 


d>Z 


(A- 10b) 


Of  course,  the  above  can  be  found  in  many  text  books  and  we  repeat  it  here  for  the  reader's 
convenience. 

Now  let  the  instantaneous  Poynting  vector  be  given  by 


S  =  E  xH  (A-l  1) 

where  all  boldface  variables  indicate  instantaneous  quantities,  if  is  a  power  density,  i.e. 
power/unit  area  and  has  the  MKS  units  of  watts/m2.  The  total  power  crossing  some  closed 

surface  at  a  given  time  is  found  by  integrating  the  normal  component  of  if  over  the  entire 
surface,  or 
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P - H  S  ds  =  ||  S •  hda 

s  s 


(A- 12) 


where  P  is  the  instantaneous  power  in  watts.  We  can  write  Is  and  $  in  the  form 
(Reference  17) 


E(jc,v,z;f)~  Rc 


E(x,y,z)e 


tax 


(A- 13a) 


H(x,y,  z;r)  =  Re 


H{x,y,7.)el(Dl 


(A- 13b) 


and  then 


E  =  E'  cos  (oi  +  £'■'  sin  ax 


(A- 14a) 


and 


— V  — >  — > 

H  =  H'  cos  cot  +  H"  sin  ax 


(A-14b) 


d  we  have  thrown  away  all  else.  Using  Equations  A-14,  the  instantaneous  Poynting 
•'tor  becomes 


-> 


S  =  E'xh'  cosL  «w+  E' x  H"+  E" x  H'  sin  an  cos  an 


+  E'xh' 


.  -> 

sin“  cot 


(A- 15) 


An  alternate  form  of  Equation  A- 15  is 
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1  (  —*  —*  — if—*  — *  —*  — >  ^ 

S  =±\E'xH'+E"xhr\+^\  £'*//*+ £*jr//'lsin2fi>r 


(A- 16) 


>  — »  — »  ^ 
E"xH"  c 


+  -[  E'xH'-  E" xH"  |cos  2(0l 
2 


or 


—4  —4  —4  — ► 

S  =  A/q+  A'i  cos  loot  +  Ni  sin  2 ox 


where 


(A- 17) 


f 


E’  x.H'+  E"xH' 


(A- 18a) 


— ►  i  /  — >  — >  — >  — ►  A 

N\  -  -j\E'x  n’-  E“x  H"  j  (A- 18b) 

(A- 18c) 

which  is  Grimes’  notation  (Reference  6). 

Once  we  time  average  Equation  A- 17  over  one  period,  T,  we  find  that  the  integrals 
over  cos  2  cot  and  sin  2(0t  are  zero.  Thus  while 


N^2 


E'xir+  E"xH' 


=  Re 


ExH 


(A- 19) 


the  terms  containing  Ni  and  N2  average  to  zero.  Then  the  radiated  power  is  also  the  time 
averaged  power,  as  usual, 
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■ad  =  l‘ms~§  MiW*) ■«*>  =  # 

A  ~  '  '  /* 

Only  the  instantaneous  total  power  keeps  the  Ni  and  N-> 

P  =  ^  S  hda 
A 

=  <ji£  Nq  hda  +  cos  2a>r^  Nyftda  +  sin  Nyhda  .  (A-21) 

AAA 


Nq- hda  .  (A-20) 


terms.  Thus 


But  neither  nor  N2,  given  in  Equations  A-I8b  and  A-18c  are  equal  to  ~  Im(  I?  x  H^), 
which  is 


(  — t 

l-Jl 

( — »  -»  ->  -- > 

1  ExH* 

E'xH"-  E" xH' 

\  J 

1  2' 

\  / 

(A-22) 


which  is  obtained  from  Maxwell  s  equations  and  their  associated  complex  conjugates. 
Thus  there  is  no  explicit  relationship  between  the  imaginary  part  of  the  complex  Poynting 
vector  and  the  instantaneous  Poynting  vector. 
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